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We study an extended Hubbard model with the nearest-neiglibor Coulomb interaction on 
the pyrochlore lattice at half filling. An interaction-driven insulating phase with nontrivial Z2 
invariants emerges at the Hartree-Fock mean-field level in the phase diagram. This topological 
insulator phase competes with other ordered states and survives in a parameter region sur- 
rounded by a semimetal, antiferromagnetic and charge ordered insulators. The symmetries of 
these phases are group-theoretically analyzed. We also show that the ferromagnetic interaction 
enhances the stability of the topological phase. 
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1. Introduction 

The discovery of the quantum Hall effect^-* and the- 
oretical approaches to the effect^' ■^^ have revealed the 
importance of topological properties of electron bands. 
In the case of the quantum Hall states, the topologi- 
cal properties of systems are characterized by an integer 
called Chern number, which is physically measured as 
the quantized Hall conductance. A non-zero Chern num- 
ber requires a breaking of the time-reversal symmetry 
typically realized by applying strong magnetic fields on 
a two-dimensional system. The quantized Hall conduc- 
tance is also topologically protected by the structures of 
the edge states. 

Recently, topological classification of the time- 
reversal-invariant band insulators in two and three spa- 
tial dimensions has been made in seminal works. ^^^^ 
The classification can be completed by using the knowl- 
edge of the single-particle wave functions alone. Two- 
dimensional band insulators are characterized by a single 
Z2 invariant, which is equal to the number of the gapless 
edge states modulo 2 protected by the time-reversal sym- 
metry. Insulators with nontrivial Z2 invariant namely Z2 
invariant = 1, are termed as topological insulators (TIs) 
or quantum spin Hall (QSH) insulators. In three dimen- 
sions, four Z2 invariants {vfu ^i^2i^z) with fa = or 1 
characterize either strong topological insulators (STIs), 
weak topological insulators (WTIs), or trivial band in- 
sulators. STIs have vq = 1 while WTIs have some of the 
{i'iU2V3) differ from and = Q. Trivial insulators are 
characterized by the index (0; 000). 

Both in two and three dimensions, the topological 
phases are typically realized in the systems with strong 
spin-orbit (SO) interaction. On the other hand, it 
has been suggested that the extended Hubbard model 
on the honeycomb lattice can generate an effective SO 
interaction from a spontaneous symmetry breaking at 
the Hartree-Fock mean-field level and results in the QSH 
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phase. Similar phenomenon has also been proposed on 
the kagomc and diamond lattices. A common prop- 
erty is that the lattice models which are semimetals in the 
single particle problem may have topologically nontriv- 
ial insulator phases caused by the Coulomb interaction. 
Therefore, these states are called as topological Mott in- 
sulators (TMIs) and we will adopt this terminology in 
this paper. We note that TMI has also been used in a 
different context. "'^''^ 

In this paper, we study the TMI on the pyrochlore 
lattice. Systems with pyrochlore lattices have been ex- 
tensively studied because they provide us with various 
quantum effects. Recently, it has been shown that a 
tight-binding model with a SO interaction and a lat- 
tice distortion generates the STI on the pyrochlore lat- 
tice. Specifically, the SO interaction yields the STI 
with (1;000) at half fining. Possibility of STI in an ac- 
tual pyrochlore iridate has also been suggested. In our 
paper, we propose that a tight-binding model on the 
pyrochlore lattice becomes STI at half filling when the 
Coulomb interaction is switched on at the Hartree-Fock 
mean-field level. In fact a STI state emerges realistically 
solely from an appropriate nearest neighbor interaction 
in contrast to other proposed TMI, where stronger next- 
nearest neighbor interactions are assumed. We mainly fo- 
cus on half filling, where the STI has a chance to emerge 
as we shall see. In addition, competitions between other 
possible symmetry broken states are important at half 
filling as well. The system is a semimetal at half and 
quarter filling without the Coulomb interaction. 

This paper is organized as follows: In §2, we intro- 
duce a microscopic tight-binding Hamiltonian on the py- 
rochlore lattice. Its ground state becomes TI when the 
SO interaction is switched on. We confirm it from the 
Z2 invariants. In §3, we consider electron correlation ef- 
fects by employing an extended Hubbard model with a 
nearest-neighbor interaction and without the SO inter- 
action. The Hartree-Fock mean-field calculation gives a 
phase diagram which contains TI in a region of the pa- 
rameter space even without the SO interaction. 
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Fig. 1. (Color online) (a) Pyrochlorc lattice with its unit trans- 
lation vectors forming a face-centered cubic{fcc) Bravais lattice. 
The unit structure is a shaded tetrahedron consiting of sites de- 
noted by 1, 2, 3,4. (b) Unit cell with graphical definitions of bij 
and dji. 



2. Microscopic Hamiltonian and Topological 
Phases without Electron Correlations 

In this section, we introduce a microscopic Hamilto- 
nian on the pyrochlore lattice that contains a topolog- 
ically nontrivial state. We consider in this section non- 
interacting systems with the SO coupling taken into ac- 
count to provide a basis of studies on electron correla- 
tion effects in the next section. We show that a STI state 
emerges solely from an appropriate nearest neighbor hop- 
ping due to SO coupling in 3D systems in contrast to the 
previous studies.^' The pyrochlorc structure is repre- 
sented by an fee Bravais lattice with a tctrahcdral unit 
cell as shown in Fig. 1(a). We choose the unit cell vectors 
to be 

ai = (2, 0, 2), 02 - (0, 2, 2) and = (2, 2, 0). (1) 
Then the reciprocal lattice vectors are given by 

bi = 62 = f(-l'l'l) and 



63 = 



(2) 



The first Brillouin zone is shown in Fig. 2(a) with high- 
symmetry lines and points. 

Our starting point is the tight-binding Hamiltonian 

Ho = -tY. (4c,. + h.c), (3) 

where cj^ creates an electron with spin a on the site i and 
(ij) denotes a nearest neighbor bond. In the momentum 
space, eq. (3) becomes 

Hq = ^ cl^,a^oik)nami3Ckm^ (4) 
knniaP 



with 
Hoik) - 

A{k) = 



-2tA{k) (g) ao, 

/O COS^hx ^y) COsi^ky - 

cos{kx - 



V 



kz) cos{kx + kz) ^ 
kz) cos{ky + kz) 
cos{kx + ky) 



/ 

(5) 

where c^„g. is the creation operator of a Bloch electron 
on the sublatticc n = (1, 2, 3, 4) with the wave vector k 



and spin a, and do is the identity matrix in spin space. 
The lower triangle of the matrix A{k) should be filled so 
that the matrix is Hermitian. The spectrum of eq. (5), 
shown in Fig. 2(b), is described analytically as 



-5(1,2) (fe) = -2t (l ± V^ + Bik) 

B{k) ~ cos 2fcj; cos 2fcj, -f cos 2^2: cos 2fc2 

-t- cos 2fcy cos 2A;2 (6) 



and 



^;(3,4)(fe) - 2t. 



(7) 



Here, E{2){k) touches the two fiat bands E(^4^^{k) 
quadratically at the F point at E/t = 2.0, while E(i)(k) 
and E^2){k) are degenerate at the diagonal line of the 
square faces of the Brillouin zone, namely along the line 
WX at E/t = -2.0. 

Next, we introduce a spin-dependent nearest-neighbor 
hopping which is crucial for the emergence of the topo- 
logical states. We consider the SO interaction of the form 



+ /i.e. 



(8) 

where A is the SO coupling strength, cr is the Pauli spin 
matrix, bij is the vector from the center of a tetrahedron 
to the midpoint of the bond (ij) constituting the tetra- 
hedron edge and djj is the vector connecting the site j 
to i. Vectors bij and dij are graphically shown in Fig. 
1(b). Since each bond participates in forming only a sin- 
gle tetrahedron, 6^ is uniquely determined when a bond 
(ij) is specified. We show that this form of the SO inter- 
action is justified from a group theoretical argument in 
§ 3.2. We note that this SO interaction may be viewed 
as a three dimentional generalization of the SO inter- 
action introduced on the kagome lattices in Refs.l4, 20, 
where the starting point of bij is the center of the tri- 
angle. However, there exists a clear distinction between 
the kagome and pyrochlore lattices: In the case of the 
kagome lattice, bij x dij is perpendicular to the kagome 
plane for all the bonds (ij). On the other hand, in our 
thrce-dimcntional system, the vectors bij x dij are not 
even coplaner, which leads to the absence of a unitary or 
an antiunitary transformation (ex. a spin rotation or a 
time- reversal transformation) that transforms A to —A. 
Therefore, the spectrum of Hq + Hso with positive A and 
negative A differs each other as shown in Figs. 2(c) and 
(d). A similar behavior was reported in the case of the 
next nearest neighbor SO interaction.^^'' In the present 
case, -ffso with positive A opens a gap at half filling while 
negative A does not. 

Now we study the topological invariants of the insulat- 
ing states derived from Hq + Hgo . We can easily evalu- 
ate (I'o; J^i 1^2^3)1 since our Hamiltonian has the inversion 
symmctry.^^^ The invariants can be calculated from the 
parity eigenvalues ^2m(ri) of the 277i-th occupied energy 
band at the 8 time-reversal-invariant momenta (TRIM) 
Ti which satisfy Ti = Ti + G, with an appropriate recip- 
rocal lattice vector G. The 8 TRIM in our system can 
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Fig. 2. (Color online) (a) First Brillouin zone of fee lattiee with 
high-symmetry lines and points, (b) Band structure of the Hamil- 
tonian Hq in eq. (3) or (5). (c) and (d) Band structure of the 
Hamiltonian Hq + H^o in eqs.(3) and (8) for A = O.lt and 
A = —O.lt, respectively 



be expressed as 

r 



i=("in2"3) 



+ 722^2 + 7l3b3)/2, 



(9) 



with Uj = 0, 1, where bi are the primitive reciprocal lat- 
tice vectors. Then is obtained by the product 



n,=0,l 



n 



(10) 



where 6i = nm=i ^2Tn(ri), with the parity eigenvalue 
^2m given below and N being the number of Kramers 
pairs of occupied bands which are paired by the time re- 
versal. If we select the site 4 in Fig. 1(a) as the center of 
the inversion then the parity operator P acts as Pckna ~ 
C-kna on the annihilation operator in the k space. Since 
C-Tiua = crino-e"*"" '"' is satisfied at the 8 TRIM with 
o.n{n = 1, 2, 3) being defined by eq. (1) and 04 = 0, par- 
ity eigenvalues for the eigenvectors of HQiTi) -\- Hso{Ti) 
are nothing but the eigenvalues of the 8x8 matrix 



Pr, = diag(e-*"i-^Se-'°2-^%e-*"-^-^Se-*°^-^^) ® ctq, 
where diag(a, 6, c, d) denotes a diagonal matrix whose 
elements are a,b,c and d. Here, _ffgo(fc) is defined so 



as to satisfy iJso = J2 
'^2m(ri) is obtained by 

PriUri2rn = ^2m(ri)ttri2m, 



nmaf) ^kna-'^Soik)naml3Ck,nf3 ■ Now 

(11) 



where Mri2m is the 2m-th eigenvector of the matrix 
Ho{Ti) + 'Hso{Ti). 

We first consider the case at half filling, where the gap 
opens for A > 0. We find that the index is (1; 000), which 
indicates that the phase is a STI for A > 0. 

On the other hand, at quarter filling, there exists a 
subtlety. We observe that Dirac type band crossings are 
present even for nonzero A irrespective of its sign at three 
X points at iJ/t = — 2. At quarter filling, the Fermi level 
is located at this Dirac point at E/t = —2. Therefore, 



even for nonzero A, the gap does not open that makes 
Z2 invariants ill defined in contrast to the case at half 
filling. To make them well defined, the degeneracy at the 
Fermi level must be removed. Therefore, at quarter fill- 
ing, each filled band must be completely separated from 
the empty bands to study physics of topological insula- 
tors. For that purpose, one may introduce a small lattice 
distortion labeled by I = 1,2,3,4 as described by the 
Hamiltonian 



Hdis ^-d-Y^ (sgn(l : ij)4^Cj^ + h.c. 



(12) 



where d is the strength of the distortion, and sgn(Z : ij) 
is + when the sublattice label of the site i or j is / and 
— , otherwise. Non-zero d opens a gap at all the Dirac 
points. Such a distortion may be realized by trigonally 
uniaxial pressure, for instance in the direction (1,1,1). 

In fact at quarter filling, d = is the critical point 
separating STI and WTI. The system is STI for positive 
d and WTI for negative d. This can be understood since 
iJo + -ffdis '^ith d = —t makes independent stackings of 
kagomc lattices. We find that {viV2V-i) depends solely on 
/ and not on the sign of d. Therefore our Hamiltonian 
can yield 8 different topological phases depending on the 
sign of d and Z = 1, 2, 3, 4. The vector G^, = 2 3 ^i^*' 

which can be defined independently of the choice of bi,^-' 
is parallel to the axis passing through the center of the 
tetrahedron and the site I in Fig. 1(a). This result is con- 
sistent with the argument in Ref.7, which suggests that 
the WTI is interpreted as a pile of two-dimensional QSH 
states stacked in the Gi, direction. The topological in- 
variants obtained here are identical to those in Ref.l8, 
when the next nearest neighbor SO interaction is consid- 
ered. Although we found that our Hamiltonian and their 
Hamiltonian are linked without closing a gap, the signif- 
icance of the present result is that the nearest-neighbor 
SO coupling in our Hamiltonian can be reproduced as 
an effective hopping matrix solely by the electron corre- 
lation of a realistic value, as we will discuss in the next 
section. 

The topological nontriviality of our system can also be 
seen by studying the surface states of the slab geometry, 
where we take the periodic boundary conditions along 
the two spatial directions and an open boundary con- 
dition along the other direction as shown in Fig. 3(a). 
Figures 3(c) and (d) show the numerical diagonalization 
of the lattice Hamiltonian Hq + Hso + -^dis '^^ ^^'^ ^^^^ 
with a positive and a negative values of d, respectively, 
where the sublattice index / is chosen in order that the 
axis passing through the center of the tetrahedron and 
the sublattice I become perpendicular to the slab. We 
plot the band energies along the lines that connect the 
four two-dimensional TRIM as shown in Fig. 3(b). The 
spectrum of the gapless surface states are clearly seen for 
both of d = 0.1 and d = —0.1, at E/t ^ 0, namely around 
the Fermi level at half filling whereas at E/t ^ —2 only 
for positive d. This confirms the above arguments. 
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(b) M2 




Fig. 3. (Color online) (a) Slab of pyrochlore lattice. We take pe- 
riodic boundary conditions along the dashed (green) lines and 
open boundary condition along the dotted (red) lines, (b) Sur- 
face Brillouin zone specified with high-symmetry points and 
lines, (c) and (d) Band structures of the lattice Hamiltonian 

Ho + Hso + J^dis "'^ ^^'^^ ^ = t'-^* = ±0.14, re- 

spectively. The sublattice index I is chosen in order that the axis 
passing through the center of the tetrahedron and the sublattice 
I become perpendicular to the slab. 



3. Interaction EfTects by Hartree-Fock Mean- 
Field Calculation 

In this section, we study electron correlation effects 
without the SO interaction. We show results of Hartree- 
Fock mean-field calculations of the extended Hubbard 
model with the obtained phase diagram under competi- 
tions among various phases. Our result shows that, even 
without an explicit SO interaction, the form identical to 
cq. (8) emerges as a consequence of a spontaneous sym- 
metry breaking in an appropriate range of the interaction 
parameters. It leads to the emergence of the STL We first 
concentrate on the half-filled case here. 

3.1 Extended Hubbard model and Hartree-Fock approx- 
imation 

The extended Hubbard model under consideration 
takes the form 

H = Ho + U'^n.^nii + V'^ninj, (13) 

where n^^ = cla^io-, = n^^ + with n^^ being the 
density operator = cl^Ci^ and (ij) is the summation 
over the nearest neighbor pairs. The onsite interaction 
is proportional to U and the nearest neighbor interac- 
tion is given by V. We decouple the on-site interaction 
according to 

«it"ii ~ n^fin^i) + - (71.;-^) (71.4) 

+ (4tC^i)(4iC.t>> (14) 
and the next neighbor interaction as 

UiUj sa n,{nj) + {ni)nj - {ni){n-j) 

- (4c,7^>(4c„)). (15) 

We note that when the system is a semimetal, where no 
symmetry is broken, {rii^) = (n^), (cj'.^Cj^) (c^j^Cj^) ^ 

0, and (c||Cj4.) = (c^Cj-i-) = are satisfied. 

3.2 Symmetry o] pyrochlore lattice 

Before listing up candidate solutions of our mean- 
field calculations, we summarize the space group of the 
pyrochlore lattice, which will be useful to identify the 
mean-field Hamiltonians of the candidate states intro- 
duced in the next subsection. The space group of the py- 
rochlore lattice is Fd3m, which contains 24 symmorphic 
elements of tetrahedral point group 43771 and 24 non- 
symmorphic elements. The non-symmorphic elements, 
which include inversion, can all be described by a com- 
bination of the translation, tetrahedral point group, and 
inversion. Therefore, we only list the 24 elements of the 
tetrahedral point group Td, by taking the position of one 
tetrahedron's center fixed. We choose this tetrahedron 
consisting of the sites 1,2,3,4 in Fig. 1(a). Its elements 
are: 

(1) 1: The identity 
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(2) C^{i = 1,2,3,4 n = 1,2): 2n7r/3 rotation around 
the axis passing through the site i and the center of 
the tetrahedron. 

(3) Ca{a, — x,y,z): tt rotation around the axis which 
passes through the center of the tetrahedron and is 
parallel to the a axis. 

(4) (Jdij{i 7^ j i,j = 1,2,3,4): Mirror reflection about 
the plane determined by the sites i, j and the center 
of the tetrahedron. 

(5) s"(a — x,y,z n = ±1): Combination of n7r/2 ro- 
tation around the axis parallel to the a axis which 
passes through the center of the tetrahedron and 
mirror reflection about the plane a = 1/2 

Then the invariant group of the TJq is described by 



Go = T X Td X {1, F} X 5C/(2) x {1, 9}, 



(16) 



where T is the translation of the fee lattice, P is the 
inversion and is the time reversal. The invariant group 
of a candidate phase must be a subgroup of Gq. 

Now we consider the invariant group of the Hamil- 
tonian when the SO interaction is taken in to account. 
Since the interaction is coupling of the orbital and spin 
angular momentum, it preserves the inversion and time- 
reversal symmetry. On the other hand, the elements of 
Td must combine with a specific element of SU{2), since 
the elements of the invariant group must rotate orbital 
and spin angular momentum cquivalently. To further dis- 
cuss the symmetry, we define the spin rotation operator 
u{n,9) which rotates the spins by around n axis. Then 
the invariant group Gso is described by 



Gso = Tx {1,P} X Lso, 
where iso contains 



(17) 



^so 



{1,©} 



1 



C^u{n,,lmT) {i = 1,2,-^,^71 = 1,2) 
X < Cau{ea,ir) [a = x,y,z) 

adiju{nki,TT) {{i,j,k,l} = {1,2,3,4}) 
.sXea,-f7r) {a = x,y,z n = ±l) 



(18) 



Here, rii is defined by 

n.^-^j(l.-l.l) 



n.. -^(-1,1,1) 



-1), n4 



1,-1,-1) 



(19) 



riij is parallel to the bond connecting the sites i and j on 
the same unit cell and is a unit vector of a direction. 
Under this symmetry, the form of the Hamiltonian is 
restricted to H = Hq + Hso defined by eqs.(3) and (8) 
as long as we consider the nearest neighbor hopping. 

3.3 Candidate phases 

Let us now introduce the phases which we find stable 
at half filling. Besides the topologically nontrivial phase, 
we also examine possible spin density waves (SDW) and 




Fig. 4. (color online), (a) Unit cell of CDW phase. Black and 
light (green) sites stand for charge-rich and charge-poor sites, 
respectively, (b) Unit cell of SDW phase. The arrows denote the 
magnetization on each site, (c) Unit cell with arrows indicating 
directions of bij X dij. These are relevant to the Fock terms of 
the SDW phase and the TMI phase. 



charge density waves (CDW). In our mean-field calcula- 
tions, we restrict to the phases which do not break the 
translational symmetry of the fee lattices. 

3.3.1 Semimetal 

For relatively weak interactions, the semimetal phase 
arising from Hq remains stable. This phase retains all 
the symmetries of the Hamiltonian Hq. 

3.3.2 Charge- density wave 

A CDW phase, where the expectation values of the 
particle number differ by sites, is expected in the limit 
of strong nearest neighbor repulsion V . Assuming that 
the full translational symmetry of the fee lattice is in- 
tact, we find that the CDW phase with two charge-rich 
sites and two charge-poor sites per unit cell is energet- 
ically favored. When we define pi with {uic) = f + pt, 
with / being the filling factor, the phase is, for instance, 
characterized by 



Pi= P2= P, P3 = Pi 



'P. 



(20) 



where p is the order parameter of the CDW phase. The 
unit cell is schematically shown in Fig. 4(a). The sym- 
metry of this state is described by 

GcDW =Tx{l,G3, adi2,^d34}x{l,P} X 5(7(2) X {1,6}. 

(21) 

This symmetry allows the bond expectation val- 
ues to have three degrees of freedom, which are 
(c1o-C2<t), (c3o-C4cr): ^nd remaining 4 equivalent bonds. 
When we define gij = 
tonian takes the form 



(Cja-Cia), the mean-field Hamil- 



H 



MF.CDW 



Hq + Y,iU - 4.V)p,n, 



- ^Yl (dijct^ja + h.C.) 

+ {-AUp^ + Wp^ 

+ AVgl^ + AVgl^ + l&Vgl:,) L\ (22) 

where is the number of the unit cells. We note that 
513 = 514 = 523 = 524 and gij = gji is satisfied under 
this symmetry. 

3.3.3 Topological Mott insulator 

TMI phase is favored at an intermediate strength of 
U and V , similarly to the QSH phases in two dimen- 
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sions^^'^^^ and the TMI phase in the diamond lattice. 
As mentioned above, the invariant subgroup of Hq + -ffso 
is given by eq. (17). The bond expectation value corre- 
sponding to this symmetry is described by 



= 9<^Q 



- iV2i 



bij X dij 
\b,j X d,j\ 



(23) 



al3 



where nonzero z reduces the invariance subgroup of the 
system from Go to Gso- The order parameter z is ex- 
plicitly written as 



a0 



bij X dij 
\bij X dij\ 



(24) 



Pa 



On the other hand, g does not break any symmetry of 
the system. The mean-field Hamiltonian for the present 
case is 



MF.TMI 



Ho 



^ C 



h.C.) 



V2Vz ^ I ic. 
{24:Vg^ + 4:8Vz'^)L^. 



bij X dij I 



■ fyot/SCjp + h.C. 



(25) 



We note that the energy remains invariant under the 
global spin rotation and directly leads to an SU{2) 
degeneracy. When z ^ 0, the invariance subgroup of 
-^MF.TMI is described by Gtmi = Gso- As mentioned 
above, the system becomes the STI for the Hamiltonian 
eq. (3) and (8). Therefore the ground state of the Hamil- 
tonian also is the STI when the solution with z 7^ is 
stabilized. 

3.3.4 Spin density wave 

When U becomes larger, antiferromagnetic phases may 
be stabilized. The phase is characterized by the order 
parameter A with, for instance, (Si) = ^tii, where 
'^i ~ '^ia'^apCip/2 is the local spin operator and rii is 
defined by eq. (19) for i = 1,2,3,4. This is nothing but 
the "all-out" phase illustrated in Fig. 4(b). This phase 
breaks the time-reversal symmetry and the spin rota- 
tional symmetry. The invariant subgroup of this phase is 
described by 



T x{l,P}x L 



SDW- 



(26) 



GsDW = 

Here isDW contains 
1 

G^u(ni,|n7r) (i = 1, 2, 3, 4 n = 1, 2) 
^SDW = < Cau{ea,TT) {a = x,y,z) 

(Jaiju{nki,TT)e {{ij,k,l} = {1,2,3,4}) 
sXea,-t7r)e (a n ±1) 

/'27) 

Using this symmetry, the form of the Fock terms is re- 
stricted to 



4^ 



= gao + 77- 



n. 



iV2z 



bij X dij 



(28) 



6.5 
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Fig. 5. Phase diagram of extended Hubbard model on pyrochlore 
lattice. The solid (dashed) lines indicate first order (continu- 
ous) transitions. SM, TMI, CDW, and SDW denote semimetal, 
topological Mott insulator, charge density wave and spin density 
wave, respectively. 



where g, rj and z are order parameters. The order pa- 
rameter z is nothing but that appeared in eq. (23) and 
this term must be included since nonzero z preserves the 
invariant subgroup Gsdw- This indicates Gsdw C Gtmi 
is satisfied. We note that nonzero A or 77 reduces the in- 
variance group of the system to Gsdw- The mean- field 
Hamiltonian of the SDW takes the form 

^fMF.SDW = Hq - C/A^cJ„(nj - <T)ai3Cip 

ia/3 



h.C.) 



+ (4t/A2 - 



rii 



■ 0'at3Cjl3 



lV2z 



I bij X dij 
\b^j X dy l 



• cra/3Cjf3 + h.C. 



2Wg'^ + 2W-q^ +ASVz^)L^. (29) 



3.4 Phase diagram 

We have searched the global minimum with respect 
to its order parameters and compared energies among 
different possible phases. The U — V phase diagram is 
shown in Fig. 5. The semimetal phase exists at a small 
U and V region. When U is dominant, the SDW is sta- 
ble and the transition from the semimetal to the SDW 
is of the first order. The CDW is stabilized when V is 
dominant and the transition from the semimetal to the 
CDW is also of the first order. Most importantly, we find 
the TMI phase at intermediate U and V , where the tran- 
sition from the semimetal is continuous. The transition 
between the semimetal and TMI occurs at V/t ~ 2.6, 
and the order parameter z grows continuously from that 
point. Other details such as critical exponents are beyond 
the scope of this paper, and would be future issues. 

The phases at other fillings are also the subject of in- 
terest. At quarter filling, we find CDW phases with three 
equivalent charge-poor sites per unit cell for relatively 



J. Phys. Soc. Jpn. 



Full Paper 



Kurita, Yamaji and Imada 7 



large V and a ferromagnet for large U limit. The system 
is a semimetal for relatively small U and V. We note that 
the topologically nontrivial states are not seen even as a 
local minimum. At three-quarter filling, the state with- 
out a symmetry breaking is unstable due to the diver- 
gence of the density of states at the Fermi level. We find 
ferromagnet and CDW with three equivalent charge-rich 
sites per unit cell when U and V dominates, respectively, 
whereas we do not find the topological insulator phase. 

3.5 Discussion 

We first discuss the possibility of experimental realiza- 
tion of the TMI along with the present mechanism. The 
calculations in §3.3 and 3.4 are for T = 0, while we can 
treat the system at nonzero temperatures as well in the 
grand canonical ensemble. For V/t ~ 3.0 we find that the 
transition to TMI occurs at Tc ~ O.Olt. 

To realize the TMI state in experiments, we propose 
that compounds like Cd20s207 and TI2RU2O7 are the 
candidates since their band structures near the Fermi 
level resemble that of our present model. Especially 
in TI2RU2O7, it is reported in Ref.23 that bands which 
consist of antibonding states of Tl 6s and O 2p orbitals 
exist near half filling. The paper also shows a calculated 
band structure of the tight-binding model for Tl 6s or- 
bitals. Though the hopping parameter of the model is 
not given in the paper, we estimate t ^ 0.3 cV from 
the band width. Then the critical temperature to TMI 
is Tc ~ 3 X 10 K, which is accessible in experiments. 
To realize the TMI state in such materials, we have to 
change the Fermi level and control the parameters in an 
extended Hubbard model. We propose that an electron 
doping as well as application of pressure are effective for 
this purpose. 

Compared to other lattice models which result in 
the topological insulators at the Hartree-Fock mean-field 
level, wc propose that the pyrochlore lattice is the best 
candidate for the experimental realization for the follow- 
ing two reasons. First, since the tight-binding model on 
the pyrochlore lattice becomes the topological insulators 
solely with the nearest neighbor hopping, we do not need 
a strong next nearest neighbor Coulomb repulsion that 
had been needed in other cases. ^'^^ Second, the tran- 
sition is possible at T 7^ 0, because our model is three 
dimensional in space in contrast to the kagomc lattice. 

It is interesting to further ask what kind of interaction 
enhances the stability of the TMI state in the compe- 
titions with the other phases. Here, we show that the 
ferromagnetic interaction of the form 

^^FM = — J ^ Si ■ Sj, (30) 

suits for this purpose. This can be easily confirmed by 
calculating coefficients of the order parameters of the 
mean-field Hamiltonians. The coefficients of the order 
parameters in eq. (25) change as 

Vg^ (v-^?jg, Vz^ (^V+^?jz, (31) 

from which we see that the interaction stabilizes the TMI 
phase since it is characterized by a nonzero z when the 



time- reversal symmetry is intact. On the other hand, the 
ferromagnetic interaction destroys the SDW state, be- 
cause it is antiferromagnetic, and does not make differ- 
ence from the CDW state, because the coefficient of p in 
eq. (22) does not change. 

In the parameter space of V and A, the Hamiltonian 
between SM and TMI is either continuous or first order. 
In the TMI phase, the spontaneous symmetry breaking 
takes place at A = 0. On the other hand, the topology 
change is also involved at the transition. The interplay 
of topology change and the symmetry breaking is an in- 
triguing future issue in terms of theory of quantum crit- 
icality.24.25) 

4. Conclusion and Summary 

In the first part of this paper, we have presented mi- 
croscopic tight-binding models on the pyrochlore lattice 
which exhibit topological insulator phases. It has turned 
out that, in tight binding models, the nearest neighbor 
transfer is sufficient to realize the topological insulator 
phase without relying on transfers between pairs at fur- 
ther distance, which is, as far as we know, the first case in 
3D systems. We have also shown that the nearest neigh- 
bor transfer to realize the topological insulator phase 
exists when SO interaction is intact. The Z2 index is 
(1; 000) at half filling, whereas at quarter filling, the sys- 
tem is on the critical point separating STI and WTI when 
distortions are absent. When we view the pyrochlore lat- 
tice as an alternative stacking of kagome and triangle 
lattice, distortions may either enhance or weaken inter- 
layer hopping transfer. According as it, the distortion 
may make the system either STI or WTI. We have also 
confirmed the topological nontriviallity of the system by 
showing the presence of the surface states on the slab 
geometry, which is consistent with the Z2 index. 

In the latter part of this paper, we have studied 
the extended Hubbard model on the pyrochlore lattice. 
We have examined candidate states at half filling from 
group theoretical arguments in the mean-field Hamiltoni- 
ans. The numerical Hartree-Fock mean-field calculations 
show that the topological insulator phase is energetically 
favored for an appropriate range of parameters with the 
nearest neighbor Coulomb repulsion. The transition from 
the semimetal to the TMI phase is continuous, while 
transitions to CDW and SDW are of the first order. 

In addition, we have also found that the ferromag- 
netic interaction enhances the stability of the topolog- 
ical state. The calculation in the grand canonical ensem- 
ble has shown that the critical temperature to the topo- 
logical phase is as much as O.Olt and is experimentally 
accessible. 
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